Finding the Kraus decomposition from a master equation and vice versa 
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For any master equation which is local in time, whether Markovian, non-Markovian, of Lindblad 
form or not, a general procedure is given for constructing the corresponding linear map from the 
initial state to the state at time t, including its Kraus-type representations. Formally, this is 
equivalent to solving the master equation. For an A'^- dimensional Hilbert space it requires (i) solving 
a first order A'^'^ x A'^'^ matrix time evolution (to obtain the completely positive map), and (ii) 
diagonalising a related x A'^^ matrix (to obtain a Kraus-type representation). Conversely, for 
a given time-dependent linear map, a necessary and sufficient condition is given for the existence 
of a corresponding master equation, where the (not necessarily unique) form of this equation is 
C , explicitly determined. It is shown that a 'best possible' master equation may always be defined, for 

approximating the evolution in the case that no exact master equation exists. Examples involving 
qubits are given. 

: 
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I ^' I. INTRODUCTION 

The theory of open quantum systems is concerned with the study of quantum systems whose dynamics is determined 
> ■ both by interactions internal to the system, and by influences that lie outside the system. These external influences 
^ ' are, of course, fundamentally quantum mechanical in nature, but under some circumstances it is useful to suppose 
^ I that they take the form of classical, stochastically varying driving terms in the Hamiltonian of the system. Either 
^H' way, what is typically dealt with is the master equation, the equation of motion for p, the density operator of the 
I— I ■ system. 

In the case in which the system of interest is interacting with an external quantum system (typically referred 
^ ' to as the environment), use of the Nakajima-Zwanzig projection operator technique [l|, Q shows that under fairly 
■ general conditions, the master equation for p — rr[x], where x is the density operator for the combined system plus 
] environment, takes the form 



m = -i [H, p{t)] + 1* jCtApi^)] ds, (1) 

o 

' where H is the system Hamiltonian in the absence of the interaction with the environment, and /Cf,s, the memory 
, kernel, is a linear map describing the effects of the environment on the system. The first term represents the coherent 
evolution of the system, while the second term gives rise to incoherent effects such as damping and fluctuations of 
the system. These effects of the environment are seen to enter via an integral over the past history of the state of the 
' system, implying that the evolution of this state at time t depends on its previous history and as such, the system 
JH is loosely said to possess non-Markovian dynamics. However, in many applications an important approximation, the 
- - Born-Markov approximation, can be made to ^ , which amounts to approximating the memory kernel above by 

^t,s[p(s)] ~ ICS{t ~ s)p{s) 

In this case, there is no memory, and the system dynamics is said to be Markovian. 

Of importance for the work reported in this paper is the fact that it has been shown via the time-convolutionless 
projection operator method [3, [J that the general master equation ([T]) can be written in a local-in-timc form, 

m^Hpit)], (2) 

where At is, amongst other requirements, a linear map such that Af(p) is Hermitean and traceless for all p. In the 
particular case of Markovian dynamics, A is independent of time (provided H is assumed to be time independent), 
and can be shown to take a specific structural form known as the Lindblad form The Markovian case is much 
studied; much less is known about the general mathematical character of non-Markovian master equations, or their 
physical interpretation. Known examples of non-Markovian master equations include those describing the damped 
harmonic oscillator and spontaneous decay of a two-level system j6|]. This is a research area that is becoming of 



2 



increasing importance with the recognition that error correction in quantum computing might not be adequately 
handled using Markovian methods, as the assumptions underlying the Markov approximation seem to be incompatible 
with the requirements assumed for error correction in quantum computing [7|. The errors that occur in quantum 
computing need not be statistically independent, i.e. the error-producing process may have memory, hence the 
need for an understanding of non-Markovicity. In [^, [lo| , non-Markovian noise models for fault-tolerant quantum 
error correction are considered. A description in terms of time-local master equations seems to be necessary for 
constructing a quantum trajectory unravelling, and perhaps a measurement interpretation, of non-Markovian master 
equations [ll [II [11, H, [Tl, [11, Ip ■ 

An im por tant feature of the derivation of the result ^ is that it relies on the existence of a certain operator inverse; 
e.g. see [l^. It is possible, at certain critical times, for this inverse to cease to exist. Some initial states evolve to 
the same final state at this time, and beyond this time there is no local-in-time master equation. This seems to 
be a common feature of the master equation description of non-Markovian dynamics. The existence of an operator 
inverse is found to reoccur in somewhat different form in the work to be described in this paper. We are able to show 
that a time-local master equation may exist even though the time evolution is not invertible, provided that certain 
consistency conditions are satisfied. 

In the case in which there is an external classical stochastic influence, the density operator is obtained by averaging 
over all the realizations of the stochastic driving term. While it is reasonable to expect that the master equation so 
obtained will be of the general form ([1]), no such proof appears to be available, though specific examples have been 
studied [18, 19, 20]. Once again this equation can be written in the time-local form 

An alternative way of representing the dynamics of an open quantum system is based on the recognition that 
the density operator p{t) must be related to its initial value p(0) by a (usually completely positive) linear map, 
p{t) = (f>t[p{0)] [2l|. This observation provides us with an alternate proof that memory-kernel master equations can 
in fact be rewritten in time-local form. In particular, consider some evolution process described by the linear map 
p{t) = (pt [p(0)] , which satisfies a memory- kernel master equation of the general form above (assuming the Hamiltonian 
evolution is absorbed into /C), i.e., 

p{t) = / ds/Cs,t[p(s)] 
Jo 

for some linear kernel map K,s,t- If it is assumed that the map 0t is invertible for the time interval considered, i.e., 
that there exists a linear map (j)^^ satisfying (f)^^ o cjjt — I where / is the identity map I{X) = X , then one can write 

p{t) = f ds {iCs.t o <i>sMo)] 

Jo 

t 

ds(/C,^to^, o^-1){^4p(0)]} 

= Mpit)], 

where At := J^dslCs,t ° 4>s ° 4>~t^^ which is of time-local form. Note that the assumption that the evolution is 
invertible is violated only if two initially distinct states evolve to the same density operator at some time t (e.g., if 
an 'equilibrium state' is reached within a finite time). This is the same kind of behaviour commented on earlier with 
regard to equation As we will show, even in this case it is sometimes possible to describe the evolution by a 
time-local master equation. Let us also note that the restriction that the time evolution be invertible need not be 
serious. Equilibrium states are typically not reached in a finite time in physical systems. Physically relevant examples 
where the time evolution is invertible for any finite time include spontaneous emission and optical phase diffusion [22| . 
The latter is also an example of a time-local master equation where the time evolution is not of standard Lindblad 
form. 

This now brings us to the issue we wish to investigate in this paper: the relationship between these two ways of 
describing open system dynamics, i.e., the relationship between A( and (pf This investigation is undertaken in two 
directions. First, if we are given a specific form for At, i.e. we are given the master equation - which need not be 
Markovian so At is not necessarily of Lindblad form - what is the corresponding linear map 4>t ? The inverse question 
is also analysed: given a linear map (pt, is there a corresponding At, and what is its form? 

We find that the first question can be answered by directly showing that given an arbitrary master equation 
p = At{p) for some linear map Af, we can explicitly construct Kraus operators {A^} 21] such that 




p{t) = = Y.ek{t)Ak{t)p{0)Al{t) 

k 



(3) 
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for a linear map (j)t, where ek{t) G {—1, 1} and where the map (f>t is completely positive if and only if ek{t) = 1 for all 
k. The construction proceeds by first selecting an orthonormal basis of Hermitian operators, which maps the master 
equation to a linear matrix equation. This is a procedure that has been used elsewhere in the analysis of properties 
of master equations [6*1 . The matrix equation is then solved to obtain (pt , and finally, a result of Choi [1^, [2^ is used 
to construct a Hermitian matrix which, when diagonaliscd (or even if only expressed as some sum of vector outer 
products) allows one to obtain the corresponding Kraus-type representation(s). An example of the method applied 
to a non-Markovian master equation is presented. 

The inverse question turns out to be more interesting. A crucial step in the analysis relies on the existence of the 
inverse of a certain matrix. If this inverse exists for all time, the construction of the master equation is straightforward, 
though not necessarily simple. However the inverse may not exist, either for all time, or at least at some critical time, 
or times. In such cases, it is nevertheless possible, provided certain consistency conditions are satisfied, to introduce a 
generalised matrix inverse, the so-called Moore-Penrose inverse, which leads to a master equation, which, however, is 
not necessarily unique. Even if the consistency conditions are not satisfied, so that no 'true' master equation exists, 
it is possible, using this procedure, to define a 'best possible' master equation. Application of the method, including 
a discussion of the possible circumstances under which the generalised approach must be followed, is illustrated by a 
simple example. 



The first step in relating a general, non-Markovian master equation p{t) = A[p{t)] to a corresponding completely 
positive map p{t) — (j)t[p{0)] relies on expressing this relationship in matrix form. The resultant equation is then the 
starting point for the two investigations outlined in the Introduction. 

Let {Ga} denote any convenient orthonormal basis set for the Hermitian operators on the Hilbert space, i.e. 



There is no need to restrict these operators to be Hermitian, but doing so gives a real representation for density 
operators. It follows that any operator X can be expressed as 



where the last equation follows via multiplication of the first by Gk and taking the trace. 

In particular, for any map ((> (though we will be concerned mainly with completely positive maps below) one has, 
expressing both p and (/"(p) as a basis sum. 



II. THE MASTER EQUATION IN MATRIX FORM 



Gi — G, 



tr[GaGb] = Sab- 





4>{p) = ^tr 4>{Y,tY[Gip\Gi) Gk ^Y.^Y[Gk4>{Gi)MGip\G, 



k, 



k I I J k,l 



which can be put in the matrix form 



0(p) = {FrfG, 



(4) 



where 



Fki:=tr[GkHGi)], n := tr[Gip]. 



(5) 



Indicating time-dependence of (f> and F via 4>t and F{t), and defining 



p{t) ■.= Mpm, 



it follows that the time evolution of p can be expressed as 



p=[F(t)r(0)pG. 



Now, suppose that p(t) satisfies a master equation of the form 



P^A{p). 



(6) 
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Note that the hnear map A may be time-dependent. To be a vahd master equation the hnear map A will have to 
satisfy certain conditions [e.g., A(p) must be Hermitian and traceless for all p], which will be assumed to be satisfied 
here (at present all one actually needs is the linearity of A). Defining the matrix L by 

Lki := tr[GfcA(Gz)] (7) 

implies immediately that this master equation can be rewritten as 

p = [Lr{t)fG. 

Note that L (and F) are necessarily real, but not symmetric in general. 
Comparing the results of the last two paragraphs gives 

Fr(0) = Lr{t) = LFr(O). 

Hence, since by linearity this equation must hold for all vectors r(0) (whether or not they correspond to a density 
operator), and noting that the G-decomposition is unique, one has the general relation 

F = LF. (8) 

This equation provides a matrix representation of the corresponding relation = A o between the linear maps cj) 
and A, and allows us to proceed in either of two directions. 

First, for a given master equation one can determine the corresponding linear map 0t. In particular, the matrix L 
is calculated via Eq. ([7]); the matrix equation ([5]) is then solved for F, subject to the initial condition F(0) = I; and 
(pt follows via the representation 

MX) = {FxfG. (9) 
In fact, solving the matrix equation formally is trivial: one has 

F(t) =F(0) rexp[ /" dsL(s)] =rexp[ /" dsL(s)] (10) 
JO Jo 

(T denotes time ordering). Note that for time- independent master equations this reduces to 

F{t) = c^* 

which can be evaluated using the standard technique of reducing L to Jordan canonical form. Details of the procedure 
to be followed to construct a Kraus-type decomposition of (f>t from knowledge of L and hence F, as per Eq. ([3|), are 
presented in section Hill 

Conversely, if we are instead given the map (f>t , we can construct the matrices F and F via Eq. ([5]) , and attempt to 
solve ([HI for L, so as to determine a corresponding master equation. This is straightforward when F is invertible, but 
leads to non-trivial considerations otherwise, and will be discussed in section IVl 



III. GENERATING KRAUS-TYPE DECOMPOSITIONS 

Choi has shown that complete positivity of (j) is equivalent to the positivity of a particular matrix S (23l . [2^ , and 
that the Kraus decompositions of (p are related to the outer product decompositions of S. The master equation gives 
us a matrix L, which in turn gives us a matrix F, characterising the linear evolution map (pt. In order to proceed from 
F to a Kraus decomposition, all we need to do is find S from F, and then diagonalise it. This procedure works just as 
well for maps which are not completely positive, leading to a Kraus-type decomposition of the form of Eq. ([3]), as will 
be shown below. In fact, our construction of S from L allows us to determine whether or not, for a given proposed 
master equation, the corresponding map is completely positive. 

The matrix S is defined as follows. Let {|q;)} denote an orthonormal basis for the Hilbert space in question, with 
a — 1, 2, . . . , iV, and for a given linear map (f> define the iV^ x iV^ matrix 

^afc := (ai|0(|a2)(/32|)|/3i). (11) 

Here the indexes a and b denote the ordered pairs (ai,a2) and (/3i,/32) respectively. We will refer to S as the Choi 
representation of (j) [23S| . 
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Noting that the Ta = |ai)(a2| form an orthonormal (non-Hermitian) basis for the operators on the Hilbert space, 
one can write (similarly as for the Ga basis) 



(b{p) = ^5ab(a2|p|/32)|ai)(/3l| 
ab 

= ^^ab|ai)(a2|p|/32)(/3i| 

ab 

= Y^^-bTaprl (12) 

ab 

It follows trivially that if (f> maps Hcrmitian operators to Hermitian operators, then S is Hermitian, i.e., 

^ba — ^ab- 

Choi further demonstrated that </> is completely positive if and only if S is positive, i.e., S > [23L [2^. 
The Hermitian property implies that one can always decompose S as a sum of outer products, i.e., 

S = Y,e,Vi^)Vi^)^ (13) 

i 

for some set of vectors For example, if the eigenvalues and normalised eigenvectors of S are given by 

SX{i) ^ \X{i), (14) 
where the are not necessarily positive as we are here not assuming (pt to be completely positive, we may then take 

Vii) = V\k\X{t), e,; = sign(A,), (15) 
corresponding to the orthogonal decomposition 

S = ^A,X(z)X(z)t 

i 

of S. Note that all decompositions of S are related by unitary maps, and hence given one, one can obtain all others 
[11,111. 

Each such decomposition of S corresponds to a Kraus-type decomposition of (f> jH,!!!. In particular, from Eqs. (fT2|l 
and (|13p one has 

0(p) = J2 ^^V{^)aTaprlV{^)l = ^e,A,;p4, (16) 

i,a,b i 

as per Eq. ([3]), where the Kraus operators are given by 

A, ^ V^(z), Ta = y(z) • T. (17) 

a 

This slightly generalises the usual Kraus representation for completely positive maps (2l| . for which case one has 
S > and hence that = 1 [ll,^. 

The above result provides us with the tool for finishing off our problem. Recall that the master equation gave us a 
matrix L, which in turn gave us a matrix F, characterising the completely positive map 4>t as per Eq. ([9]). To finally 
obtain a Kraus-type decomposition, one only needs to determine S from F, and then diagonalise it. 

Determining S from F is quite straightforward, and relies on the link between the two bases {Ga} and {tq} for 
operators on the Hilbert space. Note first that 

0t(|a2)(/32|) = $]tr[|a2)(/32|Ge] MGc) = Y.((^^\Gc\a2) MGc), 

c c 

and that 

4>fiGc) — '^^tT[GcGr]FsrGs = FgcGg. 
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Then, from the definition of S in equation (|lip . one has 

Sab = (ai|</)t(|a2)(/32|)|/3i) 

= ^(ai|F,,G,|/3i)(/32|Gc|a2) 



= $]i^.r(/?2|G,|a2)(ai|G,|/3i) 

r,s 

= J2^^rtT[GrTlG,n]. (18) 

This gives the exphcit connection between F and S as desired. 

For calculational purposes it can be convenient to consider a unitary transformation 

S(i^') := W^SW (19) 

of the Choi matrix, for some unitary matrix W. Note that this transformation defines an associated orthonormal 
basis set {Ha}, for operators on the Hilbert space, via 

Ha := ^ WabTb, 
b 

and it follows that the elements of W and S'^'^^ are given by [recalling b = (32)] 

Wab - tliHarl] = {Pl\Ha\f32), S^J^ = ^ Fsr tr [G.i? ] G.Hfc] . (20) 



The calculation of S may, on occasion, be simplified by choosing an appropriate basis {Ha}, and using equation (j20[) 
to first determine W and S^*^', to obtain the alternative decomposition 

Hp)-T.^'^rHapHl. (21) 

a, 6 

Note that the Choi matrix S is positive if and only if S^^^ is positive. 



IV. AN EXAMPLE: UNITAL QUBIT EVOLUTION 

It may be shown, noting that any orthogonal basis set for the traceless operators on a 2-d Hilbert space are related 
to the Pauli spin operators {(Jj} by some unitary transformation, that any map A( generating a time-local master 
equation for a qubit can be put in the canonical form 

3 3 

MX) = {'TfXaf -X) + Itrm ^a.af, (22) 

where :— U(i)crjU(t) for some unitary matrix U(t), and the 7j and aj are functions of time. Any Hamiltonian 
term has been absorbed by transforming to a suitable 'interaction picture' (which acts to redefine U(i) appropriately 
and leaves 7j and aj invariant). 

The evolution described by dp/dt = ^t{p) is defined to be unital if the maximally mixed state ^1 is a fixed point, 
i.e., if Af(l) = 0. Unital evolution is thus equivalent to the choice aj{t) = in the above equation. Note that, 
for unital evolution, the master equation is Lindblad-type if and only if 7j(t) > for all j (although this will not 
be assumed in what follows). With a Lindblad-type master equation, we mean a master equation that resembles 
Lindblad form, but has time-dependent decay rates Q. 

As a simple application of the methods in sections 2 and 3 we will consider the particular case of unital evolution 
where U is time- independent. By a suitable choice of qubit basis states {|1), |2)}, we can take U = 1 without any 
loss of generality, yielding the corresponding class of master equations 



p{t) = 7iaipai -I- 720-2PO-2 + 730-3Pcr3 - (7i + 72 + 72)p- 
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This class includes, for example, the physically relevant example of a two-level atom driven by a phase noisy laser, 
when no decorrelation assumptions are made |22| . In this case, the master equation is of the form 

2.9i V52 2giJ 

Here cr+ = |2)(1| and a- = |1)(2|, where |1) and |2) are the lower and upper states of the atom, and the functions gi 
and 52 are defined in terms of their Laplace transforms 

. _ s + D . {s + D){s + 4:D) + 2X^) 

~ s^ + Ds + 4A2 ' ~ s[(s + D){s + AD) + 2A2] + 2X^{s + AD) ' 

where D is a diffusion rate and A is the strength of the coupling between the two levels due to the phase-fluctuating 
classical field [l^l • 

Choosing the basis set {Ga} via Ga = Ca/V^, with (Tq := 1, the matrix L corresponding to the master equation 
follows via equation ([7]) as 



L(<) = 





-72W-73W 
-71W-73W 

-7i(t)-72(t) 



This matrix is diagonal (for non-unital maps one has Ljo — aj), and hence equation pO|) immediately yields 

F(t) 



/ 1 

Ti{t) 

T2{t) 

\o r^it) 



for the evolution matrix F, with 



where {i,j,k} is a permutation of {1,2,3}. Note that rj(0) = 1 and Tj{t) > 0. 

To find a corresponding Kraus-type decomposition, it is convenient to exploit properties of the Pauli spin matrices 
by choosing Ha = Ga in equation (|20p . so that 



r.s 



= idiag[l + Ti + r2 + Ta, 1 + Ti - - Tg, 1 - Ti + - Tg, 1 - Ti - + Tg], 



One obtains 



and hence, from equation ([21]) . that 

1 1 ^ 

Mp) - - (1 + Ti + r2 + Ta) p + - ^ [2r, + 1 - (Fi + r2 + V^)] ajpaj. 

k=l 

It follows that the evolution is completely positive if and only if the diagonal elements of S^^^ are positive. This 
is trivial for the first element, and hence complete positivity reduces to the conditions 

r, + r, <i + rfc, (23) 

where {i,j,k} run over the cyclic permutations of {1,2,3}. These inequalities, together with Tj > 0, determine a 
convex polyhedron in three dimensions, with vertices (0,0,0), (1,0,0), (0,1,0), (0,0,1) and (1,1,1). The conditions 
are satisfied automatically if "fj{t) > for all times, i.e. if the master equation is of Lindblad-type, but can also be 
satisfied in many other cases. 
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Finally, assuming that conditions (|23|) are satisfied, the density operator at time t has the explicit Kraus represen- 
tation 

3 

p(t) = ^A,p(0)4, 

with 



^0 = 




+ Ti 


+ r2- 


+ r3)^/2i 


Ai = 




+ Ti 


-r2 


-r3)i/Vi 


A2 = 




-Ti 


+ r2 


-r3)i/V2 


A3 = 




-Ti 


-r2 


+ r3)i/V3 



V. FROM LINEAR MAPS TO MASTER EQUATIONS 

We now wish to study the 'inverse' question mentioned in the Introduction, that is, given a linear evolution map 
(j)t, can we construct a corresponding master equation? This is clearly equivalent to determining whether the matrix 
equation F = LF in Eq. ^ has a suitable solution for L, for a specified matrix function F{t). In the case that F{t) 
is invertible one has the unique solution 

L = FF-\ (24) 

and hence a unique corresponding master equation. 

However, in the singular case, in which F~^ does not exist, the situation is not so clear. For example, detF = in 
such a case, and hence from Eq. ^ one requires that det F = 0. If the latter condition is not satisfied there can be 
no solutions for L, and hence no corresponding master equation. On the other hand, as we will show below, there can 
be situations for which F is not invertible for which there are a plurality of master equations which generate the same 
evolution. We will begin by discussing the general consistency conditions which must be satisfied for a time-local 
master equation to exist for a given evolution map 0t, and show how to obtain the corresponding master equation(s). 
It will also be shown that, even when the consistency condition is not met, one can still define a 'best possible' master 
equation corresponding to the map. 

A. Consistency conditions 

A time-local master equation is first-order in time. Hence, if the state is known at time t, the state at all later 
times is uniquely determined. It follows in particular that if the trajectories of two density operators p{t) and T{t) 
intersect at some time t = t' , i.e., p{t') = T{t'), then the trajectories merge into a single trajectory for all later 
times, i.e., p{t) — rit) for t > t' . This general causal requirement, that when two trajectories 'kiss' they then stay 
together forever, might be termed a 'fidelity' constraint on time-local master equations. Note that since any merging 
of trajectories is irreversible, the evolution map p{t) = (j)i[p{Q)\ is not invertible in such a case. By explicit quantum 
trajectory simulations, Diosi et al. have numerically demonstrated this kind of behaviour [25| . 

It follows that for a given map (pt to have a corresponding time-local master equation, the condition 

Mpm^M^m implies </.t[p(0)] = 0t[T(O)] for all t>t' 

must be satisfied for all density operators po a-nd tq. Noting equation Q, this may be re-expressed in terms of the 
matrix F as 

F(i') [a;(0) - 2/(0)] = implies Y{t) [x{Q) - y{Q)] = for all t > t' , 

where x{0) and y{0) denote the representations of p{0) and t(0) with respect to the basis set {Ga}- Hence, by 
linearity, any vector in the kernel of F, ker F, at some time t' must remain in the kernel for all later times. Denoting 
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the projection matrix onto the kernel of F(t) by K(i) (thus, = K and FK = 0), it further follows from the defining 
equation F = LF that FK ~ 0. Hence, the fidelity constraint reduces to the two consistency conditions 

ker F{t') C ker F{t), F{t) K{t) = (25) 

on F{t), for all times t > t' > 0. Note the first of these conditions implies that the dimension of the kernel of F 
increases monotonically, and hence provides a measure of irreversibility. Note also that if F{t) is invertible for all 
times t, then K = 0, and the conditions are trivially satisfied. In this case, the corresponding master equation is 
uniquely defined by the matrix L in equation (j24p . 

Finally, note that a third consistency condition has been implicitly assumed from the outset: that p is well defined 
at all times - otherwise, of course, no master equation of any kind can be defined. This is equivalent to the assumption 
that the matrix F{t) is well defined for alH > 0. 



B. Finding the master equation when the inverse of F does not exist 

Even if the matrix F is not invertible, so that equation (|24|) cannot be applied, a corresponding master equation 
can be shown to exist, provided that the consistency conditions ([25]) hold. 
In particular, let 

F = UDV"^ 

denote the singular value decomposition of the real matrix F. Thus, U and V are orthogonal matrices which 
diagonalise FF^ and F'I'F respectively, and D is a real diagonal matrix of the form 

D = diag[si, S2, . ■ ■ , Sr,0, . . ■ ,0], 

where r is the rank of F and si > S2 > ■ ■ ■ > Sr > 0. The singular values si, S2, . . . are the square roots of the non-zero 
eigenvalues of either of FF^ and F^F. Further, define the related diagonal matrix 

D := diag[l/si, 1/5^,0, ... ,0], 

and the projection matrix 

P := DD = DD = diag[l, . . . , 1, 0, . . . 0] = P^ 

having r Is and N'^ — r Os. The matrix 

F VDU"^ (26) 

is called the Moore-Penrose inverse of F [2^ . 

If the consistency conditions ((25|) are satisfied, a solution to the equation in Eq. ([8]), F = LF, is then given by 

L = FVDU'^ = FF (27) 

where D and F are defined above. Note that this reduces to equation (|24p when F is invertible, since F = F^^ in 
this case. 

To show there can indeed be an infinity of solutions for L (and hence an infinity of master equations corresponding 
to the evolution map (f>t), note first that VPV""^ is the projection onto the range of F, and hence 

K = I-VPV'^. (28) 

Now consider the ansatz 

L = FF + M, (29) 

for some matrix M. It follows immediately that 

LF = FFF + MF = FVDU'^UDV'^ + MF = FVPV'^ + MF = F(I - K) + MF = F + MF, 

where the last equality follows via the second consistency condition in equation ([25]) . Hence, equation ^ is satisfied 
if and only if 

MF = 0. 
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The choice M = gives the particular solution in equation ((27)) . Note that one is forced to make this choice whenever 
F is invertible. 

Finally, the explicit form of the master equation p = At(p), corresponding to equation ([27]), is prescribed by 

At(X) :=^i?,fc(i)T,XT/, (30) 

ah 

where 

Rab := ^(FF)rs tr[G,r]G,Tb], (31) 

rs 

and where the Ta were defined in section HITl To obtain these expressions, note that Rah is just the Choi matrix of the 
linear map corresponding to the solution for L, analogous to equations (fTTjl and (fT8|) . Note that one can in fact obtain 
a similar master equation for each allowable matrix M in equation Ij29p . where each such master equation generates 
precisely the same evolution map (^t (providing, of course, that all consistency conditions are satisfied). 



C. Best possible master equation when the consistency conditions are not met 

It is of interest to note that the Moore-Penrose inverse of F in equation (pS)) is always well-defined, whether or not 
the consistency conditions in section 5.1 are met. Hence, the corresponding matrix L in equation (|27[) . and the master 
equation map in (|30|) are similarly always well-defined. It turns out that, for any given evolution map 0t, the map 
A( defined by pop always yields the hest possible time-local master equation (and yields an exact master equation in 
the particular case that the consistency conditions are satisfied). 

In particular, consider the problem of attempting to model a given evolution process, p{t) :— (f>t[p{0)], by a time- 
local master equation. For example, (f)t might be defined via a particular interaction of the system with a 'bath', or 
obtained from phenomenological or theoretical considerations. This problem may be formulated more precisely as 
determining the map At for which 

— A(||= minimum, (32) 

where the norm of a given linear map (f) is identified with the Hilbert-Schmidt norm of the corresponding matrix 
representation of cf) in equation ([5|), i.e., 

m ^ ||F|| := (tr[FTF])'/' 

(note that the norm is independent of the choice of basis set {Ga})- It follows that the problem is equivalent to 
determining the matrix L such that 

||F — LFj|= minimum. 

It will be shown that the solution to the above problem, further satisfying the 'least squares' property 

||At|| = ||L|| = minimum, (33) 

is given by L = FF, corresponding to the map in equation (jSOp . This map is therefore, in this sense, the best possible 
for modelling the process 0t, even in the case where no true time- local master equation exists. 
In particular, noting that K(I — K) = and FK = 0, one has 

||F-LFjp = ||(F-LF)K + (F-LF)(I-K)||2 

= ||(F-LF)K||2 + ||(F-LF)(I-K)f 
= ||FK|p + ||F(I-K) -LF||2 
> ||FK|p. 

Further, from equations (j26p and (|28p one has FF = VPY"^ = I - K. It follows immediately that the minimum value 
is obtained for the choice L ~ FF + M, where M is any matrix which satisfies MF = 0, precisely as per equation 
(HH). Note from UPU'^F = F that M = M(I - UPU'^). Thus, noting further that FUPU'^ = F, one has 

||L|p = ||FF|p + ||M|p. 
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Hence, the 'least squares' solution is given by L = FF as claimed. 

Note that for qubits undergoing purely Hamiltonian evolution, with ihK{p) = [H, p], one finds 

n''\\Kf = iv[H\ (34) 

This suggests it may be possible to interpret the 'least squares' property (|33p more generally as some kind of 'minimum 
energy' property. 

Finally, note from section 5.2 that the norm \\(f)t — At|| vanishes identically for all times if the second consistency 
condition in equation psp is satisfied. However, this does not imply that the first consistency condition is necessarily 
satisfed at all times. In particular, the corresponding 'best possible' master equation can only track the evolution 
map (j)t from t = Q only up until times for which the first consistency condition is met: no first-order equation in time 
can model the splitting of a given trajectory into two or more trajectories. 



VI. FROM KRAUS REPRESENTATION TO MASTER EQUATION: MINIMAL DECOHERENCE OF A 

TWO-LEVEL ATOM 

Here the methods of the previous section are demonstrated for the case of a two-level atom, the evolution of which 
is described by a completely positive map of a particular form. For this example the matrix F is not always invertible. 
When it is not invertible, the consistency conditions for a master equation to exist may or may not be satisfied (and 
the master equation need not be unique). The master equation for a damped two-level atom given by Garraway (27j 
corresponds to a special case of the completely positive map we consider. We also note that this map can be obtained 
via a Jaynes-Cummings interaction with a single-mode field [285 . 

Using the notation of section 4, consider then a two- level atom (or qubit) with lower and upper levels |1) and |2) 
respectively. It will be assumed that the probability of the atom being found in the upper level at time if initially in 
the upper level at time i = 0, is given by some function p{t). Thus, p{t) is the survival probability of the excited state, 
and lies between and 1, with p{Q) — 1. No other properties of the survival probability will be assumed for now, 
other than it is differentiable. Note that exponential decay into the lower level corresponds to the choice p{t) — e~^*. 
More generally, however, p{t) may be non-monotonic, or even periodic, allowing for the possiblity of revivals. 

We will assume that p{t) is related to yo(0) by some linear map (pt- If we further allow for arbitrary decoherence 
of the off-diagonal elements of p, described by some function /(t), then linearity and the conservation of probability 
imply the form 

Hp) = P(t)p22|2)(2|+{pn + [l-p(t)]p22}|l)(l| 

+/(t)p2i|2)(l|+r(t)pi2|l)(2| (35) 

for all states p, where f{t) is some complex function satisfying /(O) — 1. The requirement that 0f maps positive 
operators to positive operators implies that det (pt (p) > for all p, yielding 

< \ f\^ det p + P22P{1 - PP22) - 1/1^22(1 - P22) 

for all states p. Taking p = |1)(1| implies |/p < p, and it is easily checked that this condition is also sufficient to 
ensure positivity. For simplicity, we will assume that strict equality holds, i.e., that f{t) satisfies 

Pit) = \f{t)\\ (36) 

This choice corresponds to maximising the off-diagonal elements of <j)t{p), i-e., to minimal decoherence between the 
levels for a given survival probability function. 

It may be checked that the above 'minimal decoherence' map 0f has the Kraus decomposition p{t) = 4>t{p) = 
j:^Ak{t)piO)Al{t) with 

Aiit) = |1 )(1|+/(0| 2)(2| 

A2it) = Vl-|/(t)P|l>(2|. (37) 



Hence the evolution is described by a completely positive map. However, it will be seen that in general the corre- 
sponding master equation need not be of Lindblad-type form. 
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A. Master equation when the inverses of (j> and F exist 



To obtain the master equation, the matrix F first must be determined, and checked as to whether it is invertible 
(and if it is not, whether the consistency conditions in section IV Al are satisfied). For simphcity, we will choose / 
to be real. It can be shown that a non-zero imaginary part of / would only add a term of Hamiltonian type to the 
resulting master equation. The basis {Ga} is again chosen as Ga — "^o-q, as per section 4. One finds 

0(Go) = V2|l)(l| + /2(t)G3, (/.(Gi) = /(OGi, 
0(G2) = f{t)G2 and ^(Gg) = /"(OGs- 

The matrix F and its time derivative F follow via Eq. ^ as 

1 

f{t) 

fit) 

Pit) -I Pit) 



fit) 

fit) 

\2fit)fit) 2fit)fit) j 



(38) 



When fit) 7^ 0, then F is invertible, with 



II Pit) 1 





l//(i) 

















(39) 



The solution for L follows via Eq. 



as 



FF"^ 







V 2/(t)//(t) 





f{t)lfit) 

fit)lfit) 

2/(t)//(t) / 



(40) 



Finally, to obtain the master equation in the form p — At(p), we calculate the Choi matrix R in equation (|3ip . 
corresponding to L, yielding 



R 



/2/(t)//(t) fit)/ fit) \ 

fit)lfit) 



V -2/(t)//(i) / 



(41) 



where the ordering {|2)(2|, |1)(1|, |2)(1|, |1)(2|} for the basis elements has been chosen. The master equation immedi- 
ately follows via equation ([50]) as 



Pit) = KMt)] = ill (2|2)(2|p(t)|2)(2| + |2)(2|p(<)|l)(l| + |l)(l|p(t)|2)(2| - 2|l)(2|p(t)|2)(l|) 
fit) 

= (2(T_p(t)a+ - a+a^pit) - pit)a+a^) . 



(42) 



It may be shown that if we do not assume that / is real, a term of Hamiltonian type is added to the master equation, 
so that it reads 



f{t)rit)+fit)f*it) 

2|/(i)P 

fit) Pit) -fit) Pit) 



(2(T-p(t)(T+ - (T+a-pit) — pit)a+(T-) + 



2|/(i)P 



[cr+CT_,p(i)]. 



(43) 
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For the case of exponential decay, i.e., p{t) = e we have f{t)/f{t) = —7/2, and then the master equation is 
of standard Lindblad form, corresponding to Markovian behaviour [6j. If f{t) is time-dependent, but in a way so 
that f{t)/f{t) is negative, then this corresponds to non-Markovian but still 'Lindblad- type' behaviour. Such a master 
equation can unravelled along the same lines as when using the usual quantum trajectory methods [ill [29l. Isol. [31I Is^ . 
as is done for instance in [l^. If f{t)/f{t) becomes positive, this corresponds to reabsorption and recoherence, and 
can be termed 'anti-Lindblad' behaviour. For this case, other quantum trajectory unravelling techniques must be used 
[IjIUjIIBI- If we have |/(ti)| < 1/(^2)1) then the time evolution map from p{ti) to ^(^2) is not completely positive. 
The total time evolution from i = to t = t2 is always completely positive, as indicated by the Kraus representation 
in (|37p . We will return to the issue of the relation between non-Lindblad master equations and complete positivity 
elsewhere. 



B. Master equation when the inverses of <f> and F do not exist 



If at some time to we have /(to) = 0, i.e., the atom occupies the lower level with probability unity, then the inverse 
of the linear map does not exist at this time. This corresponds to the inverse of the matrix F not existing. At such 
times F has the form 
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-10 0, 

implying that Fa; — for any vector with xq — 0. The projection matrix K onto the kernel of F follows as 



K(t) 



(44) 




(45) 



for f{t) — 0, and K(t) = otherwise. Hence, ker F(t) is 3-dimensional for /(t) = 0, and is 0-dimensional otherwise. 

It follows that for the first consistency condition in equation (PS)) to be satisfied, /(t) can reach zero at some time 
to only if it remains zero for all later times t > to, i.e., a time-local master equation can exist only if decay to the lower 
level is irreversible. Further, noting the form of F above, it is seen that the second consistency condition in equation 
((25)1 is satisfied at f{t) = if and only if 



FK 



/ 





fit) 








fit) 


V 






0, 



that is, if and only if / = 0, and hence, from equation ([55]) . if and only if F = 0. When these conditions are both 
satisfied, a master equation corresponding to 4>t exists even when F is not invertible. 

In particular, as demonstrated in section IVBl we may obtain a solution for L using the Moore-Penrose inverse of 
F, L = FF + M. As noted above, the consistency conditions imply that F = 0. Furthermore, the matrix F is finite 
and consequently its singular value decomposition and Moore-Penrose inverse are well-defined. Hence, the general 
solution for t > to reduces to L = M, where M is any matrix satisfying MF = 0. We therefore have 



L(t) 



fli hi ci fli 

a2 62 C2 02 

as ^3 C3 03 

a4 64 C4 04 



(46) 



where ai,bi and Ci are arbitrary time-dependent functions, for all t > to- Since (fitip) — |1)(1| whenever /(t) — 0, i.e., 
the atom is in the ground state, it follows that the master equation map corresponding to L in equation (|46p satisfies 
A(|l)(l|) = 0, irrespective of the choice of functions aj{t), bj{t) and Cj(t), with A(|2)(2|), A(|2)(l|) and A(|l)(2|) being 
arbitrary. Hence we have many possible master equations which generate the same time evolution for t > to. For 
earlier times t < to, where /(t) ^ 0, the master equation is of course uniquely specified by equation (I43p above. 

Finally, it is worth pointing out that if the 'minimal decoherence' condition (|36p is relaxed, so that the survival 
probability p(t) and decoherence factor f(t) in equation (|35p become independent (subject to the positivity requirement 
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l/P < p), then we can again calculate the corresponding master equation when the consistency conditions are satisfied. 
For the case /(io) — but p(to) 7^ 0, the off-diagonal density matrix elements vanish, corresponding to 'total 
decoherence'. In this case, the inverse of F does not exist, and the kernel of F is two-dimensional. The consistency 
conditions then imply that f{t) = for all later times, i.e., that 'decoherence is irreversible', and the set of possible 
master equations for t > to correspond to those for which A(|l)(2|) = = A(|2)(l|). If at a subsequent time ti the 
survival probability becomes equal to zero, this reduces to the case treated above, with the consistency conditions 
implying that 'decay to the lower level is irreversible'. 



C. Best possible master equation when no true master equation exists 



As shown in section 5.3, if the consistency conditions are not satisfied for the above example, then the master 
equations corresponding to L = FF -I- M, with MF = 0, minimise the quantity \\4>t — At||. Moreover, the particular 
master equation corresponding to L = FF is the 'best possible', in the sense of having the smallest possible norm. 

In order to calculate this 'least squares' solution L for the case f{t) — 0, we have to calculate the singular value 
decomposition of F in equation (|44p . as described in section 5.2. Only the first singular value is non-zero, being equal 
to v^, and the matrices U and V can be chosen as 



U = 



The Moore-Penrose inverse of F follows from equation 



F = VDU^ 



1 

V2 








V2 \ 




/I 








1 










1 











1 








1 


1 

V2 








72 / 




\0 






2B as 






/ 1/2 





-1/2 




















V 









(47) 



(48) 



and hence, using equation (|38p for F, the 'least squares' solution for L when f{t) = is given by 



/ 



FF 









(49) 



V mm -f{t)f{t) J 



At all other times, L is instead given by equation (|40p . 

If we assume that (jjt is well-defined at all times, then this implies that f{t) is finite and therefore that // = when 
/ = 0, so that L = when f = 0. In this case, for L as above, one finds that the corresponding minimum value of 
1 1 0t — At 1 1 in equation (|32p is given by zero whenever f{t) ^ 0, and by 



V2\f{t)\ 



when f{t) = 0, where here the last equality follows via equation ([55)1 . The distance between the 'true' (pt and our 
'best guess' At is therefore zero if / = whenever / = 0. This corresponds to the second consistency condition in 
equation ([25]) being satisfied at all times (for example, the case f{t) = cos^ t). 

However, as noted in section 5.3, even when the second consistency condition is satisfied identically, the correspond- 
ing 'best possible' master equation can only track the evolution map (/)f over intervals for which the first consistency 
condition in equation (j25p is also met, since no first-order equation in time can model the splitting of a trajectory 
into two or more trajectories. For example, if f{t) = cos^ t, the best possible map At cannot be used to determine 
p{t2) from p{ti) if the interval [ii,t2) includes an odd multiple of 7r/2. It is of interest to note in this regard, using 
equation (|10)) . that for t < 7r/2 the 'best possible' master equation map satisfies 



I At 



and hence diverges in the limit t 7r/2, i.e. 



'w\m/f{t)\ 

as fit) -> 0. 



2\/5| tani| 
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VII. CONCLUSIONS 



We have described a general procedure for constructing the hnear map corresponding to a given master equation, 
and vice versa. The master equation is assumed to be local in time, but may otherwise be either Markovian or 
non-Mar kovian. The method involves expressing both the master equation and the linear map in matrix form. 
Starting from a given master equation, we obtain the Kraus-type representations for the corresponding linear map 
(j}f The method works whether the linear map (j)t is completely positive or not. In fact, our construction allows us 
to determine whether or not, for a given proposed master equation, the corresponding map is completely positive. It 
will be completely positive if and only if the eigenvalues of the Choi matrix S, corresponding to the linear map, are 
nonnegative [23l . |24| . A recent paper gives a different method of obtaining the Kraus representation for the Lindblad 
equation for Markovian master equations [s^. 

Conversely, starting from a linear map 4>t, we ask whether a corresponding master equation, which is local in 
time, exists. If the (left) inverse of the map 0t exists at all times, then the construction of the master equation is 
straightforward. If the inverse does not exist at a certain time or times, then the situation is more interesting. It is 
sometimes still possible to define a master equation which faithfully describes the evolution of the system. 

If the inverse of </)( does not exist, then this corresponds to the case when some initial states evolve to the same 
state at a certain time. Thus, if one is to be able to describe the system with a first-order master equation which 
is local in time, then these states must follow the same trajectory after this point. This 'fidelity constraint' leads to 
two consistency conditions. If the map 0t satisfies these conditions, then the corresponding master equation can be 
defined, even when (j)t is not invertible. The master equation is unique if (j)t is invertible, otherwise it is not unique. If 
the fidelity constraint is not satisfied, meaning that according to the trajectories of the states which have 'merged' 
will separate again, then a true master equation, which reproduces the correct time-evolution of the system, and 
which is first-order and local in time, does not exist. However, even in this case we can define a best possible master 
equation p — At(/9), in the sense that \\(j)t — A(|| and ||At|| are minimal. As seen in the example in section 6.3, if, 
for this 'best possible' master equation, one has that ||At|| diverges at some time t', then it cannot be integrated to 
approximate (j)t after this time. 

Although master equations have previously been treated by mapping them to matrix equations Q , and the connec- 
tion between the linear map and its Choi matrix is well known [23, |2j , our intention is here to tie the two concepts 
together, and show how these tools can be used to go from a master equation to the corresponding linear map and vice 
versa. This has also led to new insights regarding the existence of time-local master equations even in the case that the 
corresponding time-evolution is not invertible. We have illustrated the methods by two examples involving two-level 
quantum systems, discussing the case of invertible and non-invertible 4>t. In these examples, we met the possibility 
of Lindblad-type behaviour, loosely speaking corresponding to time-dependent decay rates, but also 'anti-Lindblad' 
behaviour, corresponding to negative decay rates, or reabsorption and recoherence. 
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